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ION-THRUSTER PROPELLANT UTILIZATION* 
by Harold R. Kaufman 
Lewis Research Center 


SUMMARY 

The overall thruster efficiency is proportional to the propellant utilization, with 
propellant utilization defined (for singly ionized atoms) as the fraction of total propellant 
flow that is ionized and accelerated to form the ion beam. An improvement in propellant 
utilization can thus be translated directly into improved thruster performance. The 
subject of this investigation is the evaluation and understanding of maximum propellant 
utilization, with mercury used as the propellant. The primary-electron region in the 
ion chamber of a bombardment thruster is analyzed at maximum utilization. The re- 
sults of this analysis, as well as experimental data from a range of ion-chamber con- 
figurations, show a nearly constant loss for un-ionized propellant at maximum utiliza- 
tion over a wide range of total propellant flow rate. The discharge loss level of 1000 eV 
per ion was used as a definition of maximum utilization, but the exact level of this def- 
inition has no effect on the qualitative results and little effect on the quantitative results. 
There are obvious design applications for the results of this investigation, but the re- 
sults are particularly significant whenever efficient throttled operation is required. 


INTRODUCTION 

The mercury electron-bombardment ion thruster was introduced in 1960 (refs. 

1 to 2). (Although ref. 1 was published in Nov. 1960, ref. 2 was actually written first.) 
Bombardment by electrons had, of course, been used to produce ions prior to 1960. 

The need in electric propulsion was for a copious source of ions, but with a low enough 
current density to permit the integration of the source with long-lived accelerators 
operating at exhaust velocities of interest for propulsion applications. The earlier 
electron-bombardment ion sources either produced too few ions to be of interest or 

*The material presented in this report was submitted as a thesis in partial fulfill- 
ment of the requirements for the degree Doctor of Philosophy, Colorado State Univer- 
sity, Fort Collins, Colorado, June 1971. 


produced too high an ion current density to be used with practical electric -propulsion 
accelerators (refs. 3 to 10). In addition, most of the earlier sources were developed 
for specific applications, so that the source itself was often of secondary interest and 
inadequate information was given for evaluation in other applications. 

The performance of the early electron-bombardment thrusters showed advantages 
over the performance of other ion thrusters of that period (ref. 11). These advantages 
led to a substantial research program on the bombardment type of thruster (refs. 

12 to 141). Propellants other than mercury have been investigated, but only cesium ap- 
pears to be a practical alternative to mercury (refs. 142 to 156). Because the primary 
interest in space propulsion is for mercury propellant, and because there are fundamen- 
tal differences in operation with cesium, the past work with mercury will be emphasized 
in the following historical review. 


Historical Development 

The most characteristic part of an electron-bombardment thruster is the ion cham- 
ber, which started as a simple cylinder having a nearly uniform magnetic field roughly 
parallel to its axis. Most of the later improvements have been associated with the 
magnetic -field shape. First, the field windings were designed to give a decreasing 
magnetic field in the downstream direction (ref. 12). This shape is usually referred to 

as a divergent field. The strength of the optimum magnetic field was found to vary in- 

-3 -3 

versely with thruster diameter, with the optimum in the range 1x10 to 5x10 tesla 
(10 to 50 G) for a 10-centimeter -diameter thruster. An increase in magnetic-field 
strength was found to decrease ion-chamber losses at low values of field strength 
(ref. 12). The practical limit for magnetic -field strength was obtained when a further 
increase in field strength yielded a negligible decrease - or even an increase - in ion- 
chamber losses. More detailed investigations conducted later showed that ion-chamber 
noise increased with magnetic -field strength at high field -strength values (refs. 51 
and 66), which indicated that collective or wave interactions were limiting electron con- 
tainment at these higher field strengths. Iron pole pieces were added to increase the 
divergence of field lines in more recent designs (refs. 50 and 80). These pole pieces 
were originally incorporated as part of permanent-magnet designs (refs. 26 and 50) but 
have since become important tools for shaping the magnetic field. These pole pieces 
are therefore also employed on more recent thrusters that use electromagnets (refs. 

95, 112, 119, and 141). 

The foregoing constitutes the main line of magnetic -field configuration development, 
but some other concepts deserve mention. An Ioffe quadrupole magnetic field was 
adapted from fusion research (ref. 103) but gave no improvement in performance despite 
a decrease in plasma noise. The multipole ’’picket-fence” concept of plasma contain- 
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ment from fusion research has been employed in a cesium thruster (ref. 148) with a 
substantial improvement in performance. This last magnetic -field configuration has 
also been used with mercury propellant but gave no significant improvement over other 
configurations included in the same investigation (ref. 113). 

A less important, but still significant aspect of ion-chamber design is that of pro- 
pellant introduction. The first approach was to introduce the propellant at one end of 
the ion chamber and to let it flow through to the accelerator system on the other end. 

An investigation of various propellant introduction locations then showed that introduc- 
tion at the downstream end of the chamber (reverse feed) gave a substantial reduction in 
losses (ref. 37). The effects of magnetic -field shape and propellant introduction are 
interrelated, though, so that the amount of improvement possible by changing the propel- 
lant introduction mode is less for the more recent divergent-field configurations. In 
fact, for the SERT n configuration, little change was noted with introduction mode ex- 
cept for introduction near the cathode (ref. 80). Introduction near the cathode gave sub- 
stantially poorer performance. 

Considerable effort has also been expended on the development of ion-chamber 
theory. Most of the theory has been of only moderate utility, typically describing quali- 
tative trends with poor quantitative agreement (refs. 2, 14, 19, 43, 44, 51, 52, and 85). 
A thorough evaluation of published ion-chamber theory has been made recently (ref. 124). 
As concluded in this evaluation, the most detailed description of plasma processes and 
loss mechanisms has been obtained from a combination of theory and extensive probing 
of the ion-chamber plasma (refs. 21, 74, 97, and 113). To date, almost all progress 
in ion-chamber performance has been the result of cut-and-try experimentation. Ion- 
chamber theory, though, is reaching a level such that this empirical pattern may not be 
followed in the future. 

The specific area of interest for the investigation reported here is propellant utili- 
zation. Propellant utilization is defined for singly ionized atoms as the fraction of total 
propellant flow that is ionized and accelerated into the ion beam. The earliest ion cham- 
bers reached a propellant utilization of 0. 8. Shortly thereafter a scale effect was ob- 
served in which utilization was found to increase in some general manner with thruster 
size (ref. 14). Except for one attempt to correlate utilization data with a simple 
residence -time parameter (ref. 43), there were no further general advances in under- 
standing or improving utilization. Instead, the typical approach was to improve the 
utilization of a specific design at a particular operating condition by largely cut-and-try 
methods. 

The main and neutralizer cathodes have been major lifetime problems of the 
mercury-bombardment thruster. The main cathode has received the bulk of the attention 
and is considered first. The first cathode type that was used was a refractory -metal 
wire or ribbon. Ion bombardment was found to limit this type of cathode to lifetimes of 
about 1000 hours, or less, in the ion-chamber environment (ref. 29). A greater thick- 
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ness will be more resistant to sputtering, but cathode thickness is limited by conduction 
losses to supports. If the conduction losses were limited by using a heavy refractory- 
metal emitting surface heated by a fine enclosed filament, sublimation of the enclosed 
filament limits lifetime. A wide variety of matrix and dispenser cathodes (impregnated 
with strontium and/or barium oxides) were investigated next (refs. 31, 42, 53, 62, 
and 90) with the maximum lifetime exceeding 5000 hours. An oxide-coated cathode in 
which most of the surface is protected against energetic ion bombardment was also in- 
vestigated (ref. 78). A shortcoming of any cathode that depends on oxide for most or 
all of its emission, though, was brought out during the SERT n program. It is desirable 
from a reliability viewpoint to test the actual hardware to be flown before it is launched. 
A tested, but unlaunched, oxide cathode is difficult to keep from exposure to the atmos- 
phere and subsequent deactivation. This difficulty was a critical factor in the SERT II 
decision to use a hollow cathode (ref. 90). 

The mercury hollow cathode has an interesting history. A cesium plasma-bridge 
neutralizer was developed for the cesium electron-bombardment thruster (ref. 143). 

This neutralizer was a small-scale version of the autocathode used in the ionization 
chamber of the cesium thruster, and its electron emission was readily explained in 
terms of the low work function of a cesium coating. A similar design of neutralizer was 
then operated on a mercury-propellant thruster (ref. 45), where the high work function 
of mercury made the previous explanation of operation unsuitable. During the develop- 
ment of the mercury plasma-bridge neutralizer for SERT II, it became evident that the 
neutralizer cathode would easily operate at main cathode conditions. As mentioned 
earlier, the switch to the hollow cathode was then made on the basis of being able to 
ground test the actual flight hardware. A small amount of oxide is used in the hollow 
cathode, but it is buried inside the cathode and is used only to provide a few microam- 
peres of starting current. With the oxide used in this manner, repeated exposure to the 
atmosphere has no significant adverse effects (ref. 71). Another hollow-cathode design, 
now called an "enclosed keeper" design, was also considered in the early development 
of mercury hollow cathodes (ref. 45). This type of hollow cathode has recently shown 
substantial reductions over more conventional hollow cathodes in the propellant -flow 
requirements for low-emission applications (ref. 125). 

Before leaving main cathodes, one more type should be mentioned, the liquid metal 
(LM) type (refs. 56, 64, and 95). The LM cathode is similar to the hollow cathode in 
that it can be exposed repeatedly to the atmosphere when used with mercury propellant. 
The LM cathode is a workable, long-lived alternative to the hollow cathode but appears 
to require more complex solutions to the feed and isolation problems (ref. 84). It 
should also be mentioned that the baffle used to control discharge voltage in the SERT II 
thruster with a hollow cathode had its origin in a similar baffle used in ion chambers 
with LM cathodes. 
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The first neutralizers used with bombardment thrusters were simple refractory- 
metal thermionic emitters that were immersed in the ion beams (refs. 1 and 23). Such 
neutralizers are satisfactory for short tests, such as the SERT I ballistic flight (ref. 39), 
but ion erosion prevents their use for long-duration applications. Thermionic cathodes 
have also been placed near the ion beam, but such placement has always been a compro- 
mise between ion impingement from the fringe of the ion beam and high losses due to 
poor coupling to the beam (ref. 55). The plasma -bridge neutralizer avoids, this com- 
promise by providing ions to neutralize the electron space charge (refs. 83, 86, 89, 
and 131). Thus, the plasma-bridge neutralizer can be located well out of the ion beam 
and still operate with small coupling losses. 

The accelerator system that is presently used on most electron-bombardment 
thrusters consists of parallel plates with matching hexagonal arrays of holes and held 
in relative position by shadow -shielded ceramic insulators around the outside of the 
structure. Extensive experimental tests have been conducted on this type of accelerator 
system (refs. 13, 25, 41, and 61). Unlike ion-chamber theory, acceleratory theory 
has been well understood for a long time. The major cause of wear in a well-designed 
accelerator is charge-exchange erosion (ref. 22). The calculation of ion trajectories is 
well understood - although complicated and tedious. Because of the difficulty of making 
ion-trajectory calculations, early accelerator designs were created by empirical meth- 
ods. The availability of electronic computers and ion-trajectory programs, though, has 
made the theoretical approach much easier to use (ref. 105). 

Another type of accelerator system uses only one metal grid with a nonconducting 
coating on the side of the grid facing the ion chamber (refs. 54, 75, and 102). This 
composite accelerator is particularly sensitive to backsputtered metal from the ion- 
beam target and has not yet reached the reliability and durability of more conventional 
accelerators. The composite accelerator originated from discussions at the Lewis Re- 
search Center of an early version of a radiofrequency ion thruster developed by Lob 
(ref. 157). The thruster discussed had a conventional metallic accelerator, but there 
was an insulator in contact with the accelerator to take the place of the usual screen 
grid. This sharp departure from conventional electron-bombardment accelerator con- 
struction stimulated a series of experiments that ended in the construction of the first 
one-piece composite accelerators (ref. 54). 

There is a significant interaction between the ion source and the accelerator system. 
One aspect is the effect of screen-grid - open-area fraction, with the largest open-area 
fraction giving the lowest ion-chamber loss (ref. 60). Another aspect is the relation 
between ion-beam current and the current capacity of the accelerator system. For a 
given beam current and total propellant flow rate, the ion-chamber losses are lowest 
for an accelerator system operated at a small fraction of its current capacity. This last 
aspect has been shown by varying the voltage of a single accelerator configuration 
(ref. 12) and by a more general correlation of similar configurations (ref. 87). 
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Other components and aspects of electron-bombardment thrusters that have been 
investigated are feed systems (refs. 15 and 46) isolators (refs. 34, 49, and 77), arcing 
(refs. 20 and 33), and controls (refs. 28, 35, 66, and 127). The effluents of the thruster 
have also been studied, first from the diagnostic viewpoint (refs. 38 and 45) and later 
from the viewpoint of interactions with the rest of the spacecraft (refs. 72, 104, 107, 

122, 133, 135, and 140). Although the bulk of published electron-bombardment develop- 
ment has taken place in the United States, considerable work has also been done in 
Europe (refs. 63, 68, 81, 93, 98, 129, and 141). 


Present Investigation 

The objective of the present investigation is to determine the factors that control 
maximum utilization of propellant in an electron-bombardment thruster. The thrust in 
a space propulsion application is provided almost entirely from the ionized portion of 
the propellant. The un-ionized propellant thus constitutes a loss and should be mini- 
mized. The only previously published attempt to predict propellant utilization (ref. 43) 
used a residence -time approach for a neutral in the ion chamber. Although this approach 
gave the correct qualitative trends, quantitative agreement was poor. The objective of 
this investigation, then, includes obtaining a quantitative agreement better than that 
previously published. An analysis of maximum propellant utilization was published 
earlier as part of the investigation described in this report (ref. 136) and is the founda- 
tion for the theory which will be compared with experimental data. 

The thruster selected for the experimental portion of the investigation was originally 
operated with a 20-centimeter-diameter accelerator system but was masked down to a 
10-centimeter system for the work described in this report. This thruster permitted 
tests over a wide range of ion-chamber configurations, including ones with diameters 
larger than the accelerator system. An electromagnet was used instead of permanent 
magnets to permit a range of field strengths to be investigated. A strongly divergent 
magnetic field (with field-shaping pole pieces) was employed to be consistent with cur- 
rent design trends. Each ion-chamber configuration was operated over a wide range of 
propellant flow rates from near the lower limit for a sustained discharge to near the 
upper limit controlled by frequent electrical breakdowns. The thruster was operated in 
vacuum chambers located at Colorado State University and at the NASA Lewis Research 
Center. SI (mks) units are used throughout this report. 
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ION-CHAMBER THEORY 


The operation and typical values of various operating parameters are presented in 
the first part of this section. The two propellant-utilization theories that are compared 
to the experimental data are also presented, as well as a discussion of the Bohm cri- 
terion and its use in this investigation. 


Typical Operation 

A schematic of an electron -bombardment thruster is shown in figure 1. Electrons 
are emitted by the ion-chamber cathode, which is maintained 30 to 40 volts negative of 
the anode. The emitted electrons are constrained by the magnetic field (not shown in 
fig. 1) until they collide with propellant atoms. The magnitude of the cyclotron radius 
for high-velocity emitted electrons is typically 0. 1 to 0. 5 of the ion-chamber radius. 



Figure 1. - Schematic diagram of electron-bombardment thruster with power supplies. 
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Some of the ions produced by these collisions reach the right end of the ion chamber, 
where they are accelerated to form the ion beam. This ion beam is normally charge- 
density and current neutralized in space by electrons emitted by the neutralizer cathode. 
In this investigation no neutralizer was used. Instead, as is frequently done in ground 
tests, the secondary electrons from the ion-beam target were used to charge neutralize 
the ion beam. The low density and high velocity of the ions in the ion beams assure that 
only a negligible fraction of the ions recombines with electrons near the thruster. 

Un-ionized propellant atoms also escape through the accelerator and are the loss of 
particular interest in this investigation. Thruster efficiency is defined as 


Tjrjt ' 


2mP 


(1) 


where F is the thrust in newtons, m is the total propellant flow in kilograms per sec- 
ond, and P is the power in watts. In terms of the propellant flow m. that is ionized 
and ejected from the thruster and the propellant utilization m./m), equation (1) 

can be rewritten as 


TJrj. = 


vy * 2 

21 * 1 ^ 


( 2 ) 


The thruster efficiency is thus proportional to propellant utilization, and an improvement 
in propellant utilization can be translated directly into improved thruster performance. 
The understanding that might make possible such a performance improvement is the ob- 
jective of this investigation. 

During typical operation the neutral density in the ion chamber corresponds to about 
-4 18 

10 torr, or ~10 mercury atoms per cubic meter at a 500 K ion-chamber wall tern- 

17 3 

perature. The ions are only about 10 percent as dense as the neutrals (~10 /m ), and 
the electron density equals that of the ions (ref. 123). All of these densities are ap- 
proximate mean values for the ion chamber as a whole when operated at design condi- 
tions. Special locations such as the vicinity of a hollow cathode could exhibit large de- 
partures from the values given. The electron energy distribution corresponds to roughly 
a 5-eV temperature, although there are significant departures from a Maxwellian dis- 
tribution that are discussed later in this report. This approximately 5-eV temperature 
is the result of energy addition (directly or indirectly) by emitted electrons and energy 
loss by excitation and ionization. 

The Debye shielding distance for a plasma in which most of the random energy is 
in the electron population (T » T^) is 
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FF 

Z D = 7. 43X10 3 -»/ — - 
V n e 


( 3 ) 


where T is the electron temperature in eV, and n is the electron density in number 

G 0 o 

pei cubic meter. For typical ion-chamber conditions, equation (3) gives 5X10 - * 3 cen- 
timeter. This small Debye shielding distance, when compared with the ion-chamber 
diameter of 10 centimeters, justifies the assumption of n g = n^. 

Although the ion density is normally much less than the neutral density in the ion 
chamber (ref. 97), far more ions than neutrals leave through the accelerator system. 
This apparent discrepancy is due to the difference in velocities for these two species in 
the ion chamber. The neutral velocity distribution is essentially Maxwellian and is con- 
trolled by the wall temperature. The average magnitude of neutral velocity (v = 
■^8kT/7rm) is 230 meters per second for a typical wall temperature of 500 K. The loss 
rate of neutrals N q can be calculated from 


N o 


n v A, 
o o o 


(4) 


where n Q is the neutral density in the ion chamber, and A Q is the equivalent sharp- 
edged orifice area of the accelerator system. The ions, on the other hand, have a di- 
rected velocity towards the boundaries of the ion chamber, including the accelerator 
system. This directed velocity is given by Bohm's criterion for a stable sheath 
(ref. 158), 


v i 



(5) 


where q^ and m^ are the charge and mass of the ions. Expression (5) is equivalent to 
the ions having a minimum energy equal to one-half that of an electron at the most prob- 
able velocity (v m p = ^2kT/m). Expression (5) indicates a range of ion velocity is pos- 
sible. In practice, though, an equality is generally assumed. This use of an equality 
is discussed later in this section. There are corrections that can be given to the mini- 
mum velocity of expression (5). One correction that has been made is due to the depar- 
ture of electrons from a Maxwellian distribution (ref. 123). Another more difficult cor- 
rection, that apparently has been done only in a simplified form (ref. 159), includes the 
effect of ions coming from different regions at different potentials, and thus having a 
range of ion velocities. To apply the first of these corrections to the investigation of 
this report requires a completeness of probe data that has seldom been attained. To 
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apply the second would require laborious detailed calculations based on such complete 
probe data. Because of the directed nature of this Bohm velocity, the relation for loss 
rate of ions is essentially 


N i = n i v i A i 


( 6 ) 


where n. is the ion density, and is the boundary area involved in ion loss. (The 
area is usually several times larger than the neutral-loss area A Q used in eq. (4).) 
For an electron temperature of 5 eV, the Bohm velocity for singly charged mercury ions 
is 1. 55x10 meters per second. With this velocity, equation (6) indicates a loss rate of 
ions about a factor of 10 greater than that for neutrals (eq. (4)). The low fraction of 
ionization in the ion chamber is thus quite compatible with a high fraction of ions in the 
particles escaping through the accelerator system. 

The ion-chamber plasma is assumed in the foregoing to have all path lengths long 
compared to ion-chamber dimensions. The mean-free-path length for electrons collid- 
ing with neutrals or for energetic electrons colliding with slower background electrons 
is given by 


l=— (?) 

no 

where n is the density of the target species, and a is the cross section of the process. 
The path length for the total collision cross section of electrons with neutrals is about 
1 to 3 meters for an electron energy range of 5 to 40 eV. The large magnitude of this 
path length compared with ion-chamber size is the reason for using a magnetic field to 

_3 

conserve electrons. The magnetic field is weak enough (~10 T, or ~10 G), though, 

that the ion trajectories are essentially unaffected. 

The cross section for excitation or ionization by 30- to 40-eV electrons is about 
-20 

5x10 square meter (refs. 160 and 161), which gives a path length of about 20 meters. 

The electron-electron coulomb-collision cross section (90° collision) for 30- to 40-eV 
17 

electrons and a ~10 -per-cubic-meter background electron density is also about 
“20 

5x10 square meter (ref. 162), but the electron density is lower than the neutral den- 
sity so that the path length for this process is about 200 meters. Thus, the probability 
of energy loss due to excitation and ionization for a 30- to 40-eV electron is much 
greater than that due to coulomb collisions with back-ground electrons. The population 
of these energetic electrons, then, is determined by cathode emission rather than a 
Maxwellian "tail. " As a rough approximation perhaps 10 percent of the electrons are 
usually found to be energetic "primaries, " while the rest make up an essentially Max- 
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wellian distribution. * The primary electrons have been found to produce about half the 
ions in the ion chamber, while the Maxwellian "tail” produces the rest (ref. 123). This 
division of ion production is for normal operation and will be a subject of discussion for 
conditions of maximum utilization. 

Metastable states are often important in mercury discharges. Because of the low 
pressure of the discharge, though, metastable states are not important for ion produc- 
tion in a bombardment ion thruster. 


Residence-Time Approach 

The residence-time approach described in this section was published in 1965 
(ref. 43). It assumes a cylindrical ion chamber with an accelerator system covering 
one end of the cylinder and with uniform conditions throughout the chamber. The differ- 
ential equation for probability of ionization p^ is 

d Pi =v e at o^ -Pi) dn e 

where a is the ionization cross section, t is the residence time for a neutral (in the 
absence of ionization) in the volume of interest for production of beam ions, and n g and 
v g are the density and velocity of the ionizing electrons. The probability p i can be 
integrated from 0 to rjy, while the electron density is integrated from 0 to n g . This 
gives propellant utilization 77 ^ as 


Hu = 1 - ex P(- n e v e at o ) 


( 9 ) 


Thus n e v e crt 0 can be used as a correlating parameter for propellant utilization. The 
variables used in this parameter, though, are not particularly convenient. More ac- 
cessible variables will therefore be substituted. The mean electron energy can be as- 
sumed proportional to the ionization potential cp •, inasmuch as the electron energy is 
normally controlled by excitation of ground-state neutrals, and the excitation energy of 
ground-state neutrals is roughly proportional to cp Thus, the electron velocity is 

v e “ <p¥ 2 (! 0 ) 

^•This is one of the important differences between cesium and mercury electron- 
bombardment thrusters. The electron- electron coulomb -collision cross section varies 
inversely as the square of electron energy. The much lower discharge voltages used 
with cesium thus result in very rapid randomization of emitted electrons, with no dis- 
tinct primary electron group. 
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The electron density can be expressed in terms of the ion density upstream of the ac- 
celerator, which in turn can be expressed in terms of the ion-beam current density j. 
and the Bohm critical velocity. With these assumptions, one finds 


“e* 


JiW 


1/2 


<P 


1/2 


(ID 


where w is the atomic weight of the propellant. The residence time is a function of 
neutral velocity (v «: w “^/^) and ion-chamber length l , so that 

to^w 1 / 2 (12) 


From the last three relations, the correlating parameters can be written as 

n e v e CTt o° C V™ “ V ( 13 ) 

(The last step results from the propellant being restricted to mercury in this investiga- 
tion. ) Equation (9) can therefore be written as 


n v = 1 - exp(-Kj i O (14) 

where K is a constant that can be evaluated from the previously published analysis 
(ref. 43) as 5. 2. 

For equation (14) to correlate the performance of different thrusters, various con- 
ditions (such as ion-chamber wall temperature and fraction blockage of the accelerator 
system) must be nearly the same for different thrusters. In particular, the use of a 
20-centimeter chamber diameter with a 10-centimeter accelerator system is a sharp 
departure from the conditions assumed for this equation. In fact, only the 10-centimeter 
chamber data presented in this report are suitable for comparison with equation (14). 
Inasmuch as these last mentioned data will still be found to correlate poorly with equa- 
tion (14), there will be no need to modify that equation for the 20-centimeter-chamber 
data. 


Primary-Electron -Region Approach 

This analysis is directed at a condition of high propellant utilization and is con- 
cerned for the most part with primary electrons. The primary electrons do not occupy 
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the entire ion chamber, but instead are concentrated in only a portion of the volume. 

This portion is essentially restricted to the volume enclosed by the accelerator system 
and magnetic -field lines that intersect the anode. The reason for this restriction of 
primary electrons to a small volume is that the small plasma sheath on the anode does 
not reflect energetic electrons that can reach it. The existence of a primary- electron 
region is experimentally evident when detailed probe data are available (ref. 123). At 
the high-utilization condition of interest in this analysis, a high fraction of primary 
electrons is expected. The high fraction of primary electrons is in turn expected to re- 
sult in primaries causing the bulk of ionizations. The expectation of a high fraction of 
primaries is supported by probe data. 

Early bombardment-thruster designs typically had an extensive volume occupied by 
primary electrons, as indicated in figure 2. Because the region occupied by primary 
electrons in this design extends the full length of the ion chamber, a change in ion- 
chamber length will clearly affect the primary electrons. The primary-electron re- 
gion indicated in figure 2 is, of course, somewhat idealized. The finite cyclotron radius 
of primary electrons leads to a boundary with finite thickness instead of the line thick- 
ness of a limiting magnetic -field line, as indicated in figure 2. With recent divergent- 
field designs, however, the primary-electron region is localized near the screen. For 
the thruster under investigation in this report, the primary-electron region is indicated 
in figure 3. As shown in figure 3, the primary-electron region is independent of cham- 
ber geometry unless the chamber is made very short. 

The density of primary electrons increases with discharge current. For a rough 
approximation, then, the limiting condition of all electrons being primaries can be as- 
sumed for maximum utilization. Consider first the total ion production rate N from 
primary electrons. When mean values are used for the primary-electron region, 


N = n p v p n o CTV p 


(15) 


where n p and v p are the primary-electron density and velocity, n Q and a are the 
neutral density and ionization cross section, and V p is the volume of the primary- 
electron region. For the total loss rate of ions from this region, 


N = n^Ap (16) 

where n^ and v^ are the ion density and velocity toward the outer boundary of the 
primary-electron region, and A p is the area of this outer boundary. The ion velocity 
toward the outer boundary can be estimated from the Bohm criterion for a stable sheath 
given in expression (5), in which a Maxwellian electron distribution is assumed. The 
limiting condition of all electrons being primaries would probably depart from a Max- 
wellian distribution, but expression (5) should still approximate the sheath stability limit. 
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For the limiting condition of all electrons being primaries and the mean energy used for 
T , the minimum ion velocity from expression (5) is 


v i 



(17) 


As mentioned earlier, the Debye shielding distance is very small. For the primary- 
electron region, then, with all electrons in this region assumed to be primaries, 


n p = n i 


(18) 


With the substitution of equations (17) and (18) into (16) and the equating of ion loss and 
production rates (eqs. (15) and (16)), one finds 


_ V m e / 2m i 

V ^vv 


(19) 


-20 

With mercury used for im and 5x10 square meter used for the ionization cross sec- 
tion a (40-eV electrons), this becomes 


2. 34X10 16 

V /A 
P P 


( 20 ) 


A constant neutral density upstream of the accelerator system implies a constant loss of 
neutrals at maximum propellant utilization regardless of the total mass flow. 

The preceding part of this approach (eqs. (15) to (20)) has been published previously 
(ref. 136). The following extension has been made in an attempt to approximate more 
closely the effects of a non-Maxwellian electron distribution. This extension includes 
the effects of both primary and Maxwellian electrons. 

A modified Bohm criterion has been derived by Masek to give the combined effect of 
monoenergetic primary and Maxwellian electron populations (refs. 97 and 123). This 
modified Bohm criterion can be written (at the minimum ion energy) as 


v. = v 
l m 


i 


m (n„ + n ) 
e v p nr 

2m.n m 
l m 


( 21 ) 
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where v m is the most probable Maxwellian electron velocity (v m p = V^kT/m), and n 
and n m are the primary and Maxwellian electron densities. Equation (21) can be re- 
written as 



( 22 ) 


where Vp is the primary-electron velocity, and T m and & ^ are the Maxwellian elec- 
tron temperature and primary- electron energies, both in eV. Equation (18) can be re- 
written in terms of primary and Maxwellian electron populations as 


1 + n /n 

_ p m 

n. = n * 

n p /n m 


(23) 


Substitution of equations (22) and (23) into the loss-rate equation (eq. (16)) yields 


N = 

P P V 2m t f 


,T m ( 1 + n p/ n m> 


3/2 


n p/ n m 


(24) 


Equation (15) can be used again for the production -rate equation. Equating the loss and 
production rates, 


V m /2m. FT (1 + n/nj 3 / 2 

n = g- HI-™ - P --P 1 (25) 

° W ' V n m 

-20 

The mass of a mercury ion can be assumed for m^, and 5x10 square meter can be 
used for a. A value of 5 eV is reasonable for T , while ~40 eV would be expected for 
^p in a 40-volt discharge (used throughout this investigation). Equation (25) can then be 
written as 


0. 827X10 16 ^ + V n m^ 3// 

V A p V n m 


(26) 


The variation of n 0 (V p /Ap) as a function of n p/ n m was obtained from equation (26) and 
plotted in figure 4. The value of n 0 (Vp/A p ) from equation (26) is seen to be higher than 
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Figure 4. - Variation of neutral-density parameter with primary-to-Maxwellian 
density ratio - from equation (26). Primary energy, 40 eV; Maxwellian 
temperature, 5 eV. 



the same value from equation (20) for all n p/ n m < 1. Because of its inclusion of non- 
Maxwellian electron-distribution effects, equation (26) is thought to be more accurate 
than equation (20). 


Use of Bohm Criterion 


The Bohm criterion plays an important role in the preceding analysis. Some addi- 
tional comments on the use of this criterion therefore appear to be appropriate. The 
original relation given for this criterion of sheath stability was (ref. 158) 


v i 



(27) 




Figure 5. - Form of potential solutions for low-density plasma. 
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That is, the directed velocity of (cold) ions toward a collisionless sheath must be equal 
to, or greater than, the critical value of V^T e / m i' * n ma king ion-flux calculations 
for a sheath, however, it is customary to assume a simple equality. Perhaps the best 
justification for such an assumption is that it gives results in agreement with experiment 
(refs. 97 and 123). Another assumption in the preceding analysis is the use of the Bohm 
criterion at the upstream boundary of the primary-electron region. This boundary is 
not a sheath in the ordinary sense. It is, instead, a boundary that separates the 
primary-electron region with its high ion production rate from the rest of the ion cham- 
ber with a much lower ion production rate. This assumed difference in ion production 
rate is supported in the RESULTS AND DISCUSSION section by Langmuir probe data. 

To see the effect of localizing the ion production region, consider figure 5. In fig- 
ure 5(a) the ion production extends throughout the region between two walls and shows 
the customary gradual transition from plasma to sheath. In figure 5(b) the ion produc- 
tion has been localized to a central portion of this region. The solution shown in fig- 
ure 5(b) has been obtained analytically (ref. 163). The ions were found in this analysis 
to have close to the critical velocity of VqiT e / m i throughout the two regions of nearly 
uniform potential shown in figure 5(b). From figure 5(b), then, one would expect the 
Bohm criterion to apply to the boundary of the primary-electron region as long as the 
bulk of ion production takes place in this region. 


APPARATUS 

The mercury propellant thruster used in this investigation is shown in figure 6. The 
magnetic field is determined by a ring-shaped soft-iron pole piece near the accelerator 
and a cylindrical soft-iron pole piece near the cathode. These pole pieces are connected 
by four soft-iron paths (each about 3. 2 by 12. 7 mm, in cross section) on which the field 
windings are mounted. The cylindrical pole piece has about the same cross section as 
the total of the four paths and also has a field winding for most of its length. There 
were about 5. 3 turns per centimeter for all of the field windings. The construction of 
the four soft-iron paths permitted axial motion of the cylindrical pole piece for initial 
optimization. Thereafter, the cylindrical pole piece was held fixed. The data presented 
in this report were obtained with the central pole piece 3. 5 centimeter from the screen 
grid. This axial location gave good performance over a wide range of utilization. Ex- 
perimental data were taken at various magnetic -field currents, but most of the data 
were obtained at a current of 1 ampere. At this current the typical magnetic field in the 

_Q 

ion chamber (assumed to be midway between probes A and B in fig. 7) was about 3x10 

Q 

tesla (30 G). This field was somewhat nonlinear with current and rose to about 5x10" 

“3 x 

and 9x10 tesla for currents of 2 and 5 amperes. 
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Figure 6. - Thruster with adjustable ion chamber. 


CS- 58226 


A screen with 71 percent open area was used to reduce ion recombination on the 
screen. Only holes completely within a 10-centimeter -diameter circle were used to 
accelerate ions, and there were 429 such holes in the accelerator design used. The 
screen-accelerator spacing was 2 to 2. 5 millimeters over the beam area. A tungsten 
wire cathode was used to avoid both (1) the uncertainty in propellant utilization that can 
occur when an oxide cathode is inadvertently overheated, and (2) the effect of localized 
propellant concentration associated with a hollow cathode. 

Ion-chamber probes were made by using 0. 75-millimeter tantalum wire inside 
aluminum oxide tubes with an outside diameter of 1. 6 millimeters. The locations of 
the probes are shown in figure 7. The two probes closest to the accelerator (A and B 
in fig. 7) were constructed with 0. 75 millimeter of tantalum wire extending beyond the 
aluminum oxide. Probe C, because of the lower plasma density in that location, was 
constructed with a 1. 5-millimeter extension. Probe location selections were made 
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Figure 7. - Probe locations in ion chamber. 


based on published probe surveys (refs. 85 and 97). Probe A intersected a magnetic- 

field line near the edge of the cathode and was thought to be representative of the 

primary-electron region. Probe B intersected a magnetic -field line at the boundary of 

the primary-electron region, while probe C was well removed from that region. Quartz 

tubes were used for one set of probe data in an attempt to extend the range of the probe 

data to higher propellant mass flows by reducing outgassing. Tantalum wires of 0. 5- 

millimeter diameter (in holes of approximately the same diameter) were mounted flush 

with the ends of the quartz tubes in this set. 

The thruster was first operated in a 1. 2-meter-diameter, 4. 8-meter-long Space 

Propulsion Research Facility located at Colorado State University. Most of the data, 

including all the probe data, were obtained during later operation in a 1. 5-meter- 

diameter, 4. 2 -meter-long vacuum chamber at the NASA Lewis Research Center. 

Sketches of these two facilities are shown in figure 8. Both facilities used liquid- 

— 

nitrogen-cooled liners and oil-diffusion pumps to attain pressures in the 10 torr range. 
In both cases the thruster was mounted in the bell jar shown at the left end of the facility. 
Existing power supplies were used for the various thruster functions at both facilities, 
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and the propellant flow at both locations was measured by timing the drop of a mercury 
column in a calibrated glass tube. 

The meters used in voltage and current measurements are accurate within per- 
cent of full scale. Typical part-scale measurements should therefore be within ±5 per- 
cent. Because of the electrical transients that occur in other circuits when the acceler- 
ator system arcs, frequent calibration is necessary. The glass propellant-flow tube 
has a bore variation of 2 to 3 percent. This variation, together with a maximum of about 
1 percent uncertainty in column-height measurement, determines the accuracy of the 
mercury-flow measurement. 


PROCEDURE 

The thruster used was unusual in that a thruster body suitable for a 20-centimeter- 
diameter accelerator system was used with only a 10-centimeter system. Because of 
this thruster design, the ratio of thruster mass to discharge power was large. This 
large mass resulted in the slow starting transient shown in figure 9. The 0. 5-ampere 
discharge was started at time zero in figure 9. The high voltages were applied 2 min- 
utes after the discharge was started, but did not reach the values of 3000 and -1500 
volts until 12 minutes after the discharge was started. Condensation of propellant was 
assumed to have taken place on internal thruster surfaces before the thruster discharge 
was started. Evaporation after the ion-chamber discharge was started would then ac- 
count for the rapid rise and large peak in ion-beam current. Equilibrium operation was 
attained when the condensed propellant was all evaporated, which apparently took about 
1^ hours. The bumps in the descending curve of figure 9 were typical and apparently 
resulted from different parts of the thruster becoming warm enough to vaporize the con- 
densed propellant at different times. 

There was also a short-term transient effect whenever discharge conditions were 
changed. This effect was found after the initial warmup and was presumably due to the 
large ratio of internal surface area to beam area. Because of the adsorption of propel- 
lant on these surfaces at pressures of the order of 10"^ torr, a few minutes were re- 
quired to reach equilibrium when the ion-beam current was changed. The associated 
transient effect is shown in figure 10, where a hysteresis loop is shown that was obtained 
by rapidly changing the discharge current (one data point per min). Some data in the 
first paper of this investigation (ref. 136) were taken without sufficient operation at each 
condition to avoid this effect. For the data presented here the thruster was operated 
for at least 5 minutes at any conditions that involved a substantial change in beam cur- 
rent from the previous operation point. 

A transient thermal effect was also evident in the propellant flow measurement. A 
mass flow was set by holding the vaporizer constant (±1° C) at some temperature. The 
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flow was then determined by measuring the drop in height of a mercury column and the 
time required for this drop. To assure an accurate measurement of this column drop 
despite meniscus effects in an approximately 1-millimeter bore, a drop of about 4 cen- 
timeter or more was used at any vaporizer temperature. To obtain reproducible read- 
ings after any change in vaporizer temperature, though, a wait of about 15 minutes was 
found necessary. 

All steady- state operation was at an anode potential of 3000 volts and an accelerator 

grid potential of -1500 volts, with both measured relative to the vacuum facility. The 

acceleration potential difference of 4500 volts was sufficient to handle ion-beam currents 

o 

of at least 150 milliamperes without any substantial rise in impingement. No neutrali- 
zer was used, and charge -density neutralization of the ion beam was obtained by second- 
ary electrons from ions striking the various portions of the facility. 

The discharge chamber was operated at a potential difference of 40 volts for all data 
presented herein. A range of cathode emissions was obtained for each of several mass 
flow rates. The cathode emission was limited on the high end of the range by extinguish- 
ing the discharge, reaching the space-charge emission limit of the cathode, or simply 
approaching the heater power limit of the cathode. The space-charge limit on emission 
is reached when an increase in heater power results in little or no increase in emission. 
A check point of one of the initial points was usually taken after covering a range of 
emission currents. This procedure provided a spot check on constancy of other vari- 
ables while emission-range data were obtained. 

Propellant flows were limited on the low end by inability to sustain a discharge at 
40 volts and at the high end by repeated electrical breakdowns between the screen and 
accelerator grids. This range of propellant flow (and emission) was investigated with 
each of the ion-chamber configurations. 

The ion-chamber configurations were limited to three different lengths (7. 2, 10. 8, 
and 14. 4 cm) with a 10-centimeter diameter and two different lengths (9. 4 and 14. 4 cm) 
with the larger 20-centimeter diameter. These configurations provided a range of about 
4:1 for ion-chamber wall area and about 8:1 for volume. The propellant introduction to 
the ion chamber was controlled by baffles of polyimide sheet. Most of the data were ob- 
tained in the usual mode of straight-through flow from one end of the ion chamber to the 
accelerator on the other end. This mode is called "forward feed. " For comparison, 


^The current to the accelerator is usually called "impingement" in mercury 
thruster work, because a low secondary- electron emission (10 percent or less) is the 
only difference between the charge carried by ions striking the accelerator and the 
measured current. In cesium thruster work it is customary to refer to this current as a 
"drain. " This is because a cesium coating on a hot accelerator electrode can result in 
electron emission currents much greater than the ion impingement current. 
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one configuration was operated with ’’reverse feed, ” where the propellant was introduced 
through the gap between the ring-shaped pole piece and the anode. 

Most of the data were obtained at a magnetic -field current of 1 ampere, which pro- 
vided a compromise between good performance at high mass flow and ability to operate 
at low mass flow. Some data were also obtained at higher fields, though, for compar- 
ison. 

Ion-chamber probe runs were made after runs of the same configuration without 
probes. This procedure permitted the time-consuming probe data to be obtained at a 
few preselected operating conditions, instead of determining the range of operation at 
the same time as obtaining the probe data. It also permitted a check on overheating of 
the probe insulation through the effect of the resultant outgassing on beam current and 
comparison with ion-chamber performance obtained without probes. Early runs showed 
that the probes became coated with a conducting layer after 2 to 3 days of operation, 
which led to erroneous probe data. To avoid this probem, probe insulators were dis- 
carded after 1 day of operation (in addition to being checked with an ohmmeter after each 
run to verify that there was no conducting layer). The use of aluminum oxide tubes led 
to a propellant-flow limitation where (at high discharge power) outgassing occurred be- 
cause of overheating. In an attempt to extend the probe data to higher flows, the alum- 
inum oxide tubes were replaced with quartz tubes. Although this substitution avoided the 
outgassing problem, the higher cathode sputtering at higher flows led to more rapid 
coating of the probe insulators with conducting layers. Thus, no significant extension 
of the probe data was obtained by this substitution. This probe limitation is apparently 
avoided with oxide-coated and hollow cathodes, inasmuch as probe data have been ob- 
tained at higher neutral flows with these cathodes (refs. 21, 74, 97, and 113). 

In addition to the data checks described, the effects of multiple ionization and neutral 

backflow from the vacuum facility were also considered, though not included. A study 

of multiple ionization has shown that the ion-chamber discharge voltage and the cross 

sections for single and double ionization can be used to predict the fraction of doubly 

ionized atoms (ref. 16). The cross section for single ionization by a 40-eV electron is 
-20 

about 5x10 square meter, while the cross section for double ionization by an electron 

-21 

of that energy is about 1x10 square meter (refs. 160 and 161). The fraction of doubly 
ionized atoms should therefore be about 0. 02. The appearance potential for triply ion- 
ized mercury is over 63 eV; hence, no triply ionized atoms would be expected. The 
reasons for this experimental result are that most beam ions come from collisions of 
primary electrons, and that these ions leave the ion chamber so rapidly after being ion- 
ized that there is little probability of additional ionizations taking place while they are 
leaving. Hollow cathode thrusters, which have a substantial potential difference across 
the baffle system, would be expected to have a smaller fraction of doubly ionized atoms 
for the same discharge voltage. 
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In order to calculate the backflow of neutrals into the ion chamber from the vacuum 

-fi 

facility, a typical ion-gage reading of approximately 3x10 torr can be used. Most of 

the indicated reading should be due to propellant atoms during thruster operation, and 

the calibration factor (gage reading divided by actual pressure) is about 3 for mercury. 

16 

A particle density of approximately 2x10 per cubic meter would therefore be expected. 
A temperature of 300 K should be a conservative mean for the liquid-nitrogen- cooled 
tubes, the warmer baffles extending from these tubes, the room temperature end walls 
of the tank, and an ion-beam target that is warmer than room temperature. This density 
and temperature, together with an estimated effective backflow orifice area of 0. 00196 

3 

square meter, permit calculation of the backflow. This backflow, if singly ionized, is 
equivalent to only approximately 0. 3 milliampere. 

The theory of plasma probes and their use in ion thrusters has been adequately de- 
scribed in literature (refs. 21 and 41) and is repeated here only in summary form. The 
basic assumption generally used in analyzing ion-chamber probe data is that the electron 
distribution is composed of monoenergetic primaries superimposed on a Maxwellian 
background population. This assumption leads to primary and Maxwellian probe cur- 
rents of 


J 


P 



(1 + V/^ p ) 
4 


(28) 


J m = n m q A VqT m /27rm e exp(V/T m ) 


(29) 


where n and n are the densities of primary and Maxwellian electrons, q and m 

p m c 

are the charge and mass of an electron, g and T m are the energy of the primary 
electrons and temperature of Maxwellian electrons in eV, and V is the probe potential 
relative to the plasma potential in volts. The range of interest for the probe potential is 
- g < V < 0. For the usual case of S p » T m , a probe potential V approaching - <£ p 
will yield only a linear variation of primary -electron current with probe potential 
(eq. (28)). This region of nearly linear variation is shown in the plot of sample probe 
data in figure 11(a) and is approximated by the dotted line. Note that the linear region 

3 An unpublished analysis by W. R. Kerslake indicated an effective flow factor of 
about 0. 5 for the open area of a SERT II accelerator grid. The thruster used in this 
study had the same screen and accelerator grid thicknesses and hole diameters and thus 
should have the same flow factor. The accelerator open area was 50 percent of a 
10-cm -diameter beam, so that this factor of 0. 5 resulted in an effective sharp-edged 
orifice area of 0. 00196 m^. 
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Figure 11. - Typical Langmuir probe characteristics. 
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is actually somewhat curved. This curvature is due to the non-Maxwellian electrons 
actually having an energy distribution instead of being monoenergetic. The assumption 
of monoenergetic primaries, though, greatly simplifies the data analysis without de- 
parting from a reasonable approximation of the actual probe data. The contribution of 
the primaries (as approximated by the dotted line) can be subtracted from the measured 
probe currents to give a corrected current. The uncorrected and corrected probe cur- 
rents are both shown in the semi logarithmic plot of figure 11(b). It is evident that the 
lower portion of the corrected curve closely approximates the straight line expected for 
a Maxwellian distribution. The upper portion of the corrected curve is ideally a hori- 
zontal line that indicates the plasma density, and the intersection of the two straight-line 
segments indicates the plasma potential. In practice, though, a probe potential above 
plasma potential leads to an enlarged plasma sheath around the probe, so that large pos- 
itive voltages actually extract electrons from an area larger than the probe surface 
area. The plasma potential is ideally formed at the intersection of two precise straight- 
line segments. There is usually some noise in the ion chamber, and the superposition 
of curves that results while taking steady-state data is the reason customarily given for 
the rounding that occurs in this region (ref. 164). 

RESULTS AND DISCUSSION 

The ion-chamber performance data that were obtained without probes are presented 
first. Then, data for maximum propellant utilization are plotted and compared to the 
two theories presented earlier. Next, the ion-chamber probe data are presented and 
used to make a more detailed evaluation of the primary-electron-region approach to 
maximum propellant utilization. Last, a comparison is made with SERT-IT throttling 
data. 


Ion -Chamber Performance 

The performance data for five basic ion-chamber configurations are shown in fig- 
ure 12. Most of the data were obtained with a magnetic -field current of 1 ampere, as 
indicated earlier. Some of the data, though, were obtained at higher magnetic -field 
currents and are also included in figure 12. In addition, data for one configuration are 
included for operation with reverse feed, as opposed to the usual forward feed. 

The shapes of the performance curves shown in figure 12 are typical of electron- 
bombardment thrusters. Each curve has a relatively flat portion below the "knee. ” 

The location of the knee and the nearly vertical portion above the knee depend on propel- 
lant flow rate. The curves of figure 12 can be used to show the relative importance of 
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(a) Forward feed, 1 -A magnet, 10-cm diameter, 7.2-cm length. 
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Propellant utilization, 

(b) Forward feed, 1-Amagnet, 10-cm diameter, 10.8-cm length. 
Figure 12. - Ion-chamber performance. 
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(g) Forward feed, 2 -A magnet, 10-cm diameter, 7.2-cm length. 
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(h) Forward feed, 2-A magnet, 20-cm diameter, 9.4-cm length 
Figure 12. - Continued. 
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the various parameters that were varied. Comparisons of performance for the nearly 
vertical portions of the curves are adequately covered by subsequent parts of this sec- 
tion. Some comparison of the nearly horizontal portions below the ’’knee” is also of 
interest, however, and is covered here. Data for the five different configurations are 
plotted in figure 13. The neutral flows for these data are equivalent to a current range 
from 142 to 166 milliamperes (the current that would be measured if each propellant 
atom carried one electronic charge). This small variation has little effect on the per- 
formance below the knee, so that performance in this region can be compared directly. 
The variation between configurations is small below the knee - only some ±15 percent. 

In addition, there is no clear trend with ion-chamber volume or surface area. One can 
conclude that variations in ion-chamber configuration outside of the primary-electron 
region have no substantial effect on performancebelow the knee. Comments by Poeschel, 
Ward, and Knauer (ref. 113) indicate that knowledgeable workers in the ion-thruster field 
would expect such a result. Also, a ±15 percent variation can be sufficient to meet, or 
not meet, specific performance goals, and thus may appear quite important. 

One might expect a discharge-loss sensitivity to changes in the primary -electron 
region. This expectation is borne out by the lower losses found at higher magnetic -field 
currents. The higher magnetic field decreases the losses, of course, by making it more 
difficult for the high-velocity primaries to reach the anode. The use of reverse feed 
gives a small further decrease in losses, presumably through a change towards a slightly 
more uniform distribution of ion production in the primary-electron region. 

Double-valued ion-chamber performance data were not found, except possibly for a 
small region at high utilization with high magnetic -field currents. An intermittent short 
was found in the magnetic -field winding, and the extended double -valued performance 
data described in an earlier publication (ref. 136) were not encountered after this short 
was fixed. 


Maximum Propellant Utilization 

The widest range of operation was obtained at a magnetic -field current of 1 ampere. 
For this field current, though, no clear-cut maximum utilization was found (see figs. 
12(a) to (f)). Instead, the utilization above the knee increased slowly with increasing 
losses. As a somewhat arbitrary choice, the utilization at an ion-chamber loss of 
1000 eV per ion was used as a maximum. The use of some other high level of loss would 
have had no effect on the qualitative results presented and little effect on the quantitative 
results. Operation at higher magnetic -field strength gave nearly vertical lines for oper- 
ation above the knee, so a similar arbitrary choice was not needed. For consistency, 
though, the same 1000 -eV-per -ion definition was used. 
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Figure 13. - Comparison of performance for different ion-chamber 
configurations. Forward feed, 1-A magnet. 
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Residence-time parameter, Jjl,A-m 


Figure 14. - Maximum-utilization data for lU-cm-diameter ion chambers 
plotted against residence-time parameter. Forward feed, 1-A magnet. 
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Having arrived at a definition of maximum utilization, it is now of interest to com- 
pare maximum-utilization data with the predictions of the previous analyses. Maximum- 
utilization data are plotted in figure 14 against the product J^Z (ion-beam current times 
chamber length) for the three forward-feed 10-centimeter-diameter configurations at 
1 ampere magnetic field. The product JLZ is proportional to the residence-time cor- 
relation parameter Kj^Z in equation (14). The data for the three configurations were 
not correlated by the parameter J.Z, inasmuch as the short-chamber data appear to be 
shifted to the left, and those of the long chamber shifted to the right. Indeed, for the 
data in figure 14, it is clear that J- alone would have been a better correlation param- 
eter than J-Z . This result should be expected for the highly divergent -field design used 
in these three configurations, where the primary-electron region is located close to the 
accelerator system. Older thruster designs, in which the primary-electron region ex- 
tended the length of the ion chamber, were more suited to the residence-time approach. 

Maximum -utilization data are plotted against total propellant mass flow in fig- 
ure 15(a) for five forward-feed configurations operated at a magnetic -field current of 
1 ampere. The un-ionized propellant mass flow was calculated for the data in fig- 
ure 15(a) and plotted in figure 15(b). This calculation is defined by 


'°lu n -ionlzed“ (1 '” U " max> ‘'o 


total 


(30) 


From figure 15(b), the maximum -utilization data correspond to a constant un-ionized 
loss rate equivalent to 52 milliamperes, ±8 percent. This result has considerably less 
scatter than the previously published maximum -utilization data (ref. 136). The reduced 
scatter is felt to be due to improved experimental techniques, particularly those tech- 
niques associated with the long time constants of the thruster used. 

A difference between figure 15 and the results of the earlier portion of this investi- 
gation (ref. 136) is the agreement between data with 10- and 20-centimeter chamber 
diameters. There was poor agreement between data for these two diameters in the ear- 
lier results, and the lack of agreement was felt to be due to poor similarity for the re- 
gion where primary electrons reach the anode. The short 10-centimeter-diameter sec- 
tion of the anode (in the 20-cm-diameter chamber) was only 0. 6 centimeter long for the 
earlier data of the large chamber. This length was increased to 2. 5 centimeters for 
the data presented in this report (see fig. 5) to assure a dimension greater than a 
primary-electron cyclotron radius. The close agreement of data for the two chamber 
diameters in figure 15 indicates the reasoning given earlier (ref. 136) for the poor cor- 
relation was probably correct. 

The effects of a higher magnetic field and reverse feed of propellant are shown in 
the maximum -utilization data of figure 16. Reverse feed gave slightly higher maximum 
utilization of propellant, but again indicated a nearly constant loss of un-ionized propel- 
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(a) Maximum propellant utilization. 



Propellant flow rate, mA equiv 


(b) Un-ionized propellant loss. 


Figure 15. - Correlation of maximum utilization with total propellant 
flow rate. Forward feed, 1-A magnet. 


lant, but again indicated a nearly constant loss of un-ionized propellant (about 6 to 14 
percent lower than fig. 15). The data for the 2-ampere field current (forward feed) 
agreed more closely with the data shown in figure 15, but also showed a slightly lower 
loss of un-ionized propellant (0 to 10 percent lower). 

The data for a magnet current of 5 amperes in figure 12 were not included in fig- 
ure 15. At this high magnet current the discharge generally extinguished before a high 
emission could be reached. Also, a high impingement current was obtained well below 
the approximately 150-milliampere limit for lower magnet currents. This latter condi- 
tion indicated a beam-current profile that was more peaked, presumably because the 
discharge was concentrated nearer the ion-chamber axis by the high magnetic field. 
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(b) Un-ionized propellant loss. 


Figure 16. - Effect of feed and magnet current on correlation. 
Symbols same as figure IS, except for additions noted. 


As a general conclusion, then, a highly divergent-field thruster shows a nearly con- 
stant loss of un-ionized propellant at maximum utilization. This result can be obtained 
from either equation (20) or equation (26) with the proper selection of constants. In 
agreement with either of these equations of the primary-electron-region approach, the 
changes outside of the primary -electron region have no significant effect on the loss of 
un-ionized propellant. Further, even changes that do affect the primary -electron region 
can have a small effect on this value of propellant loss if the magnetic -field shape is con- 
stant (which is substantially true of the 1- to 2-ampere range included in figs. 15 and 16). 
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Ion -Chamber Probe Data 


Thorough Langmuir-probe surveys have been made of current divergent -field ion- 
chamber designs at certain selected operating conditions (refs. 97 and 113). The objec- 
tive in obtaining the probe data of this investigation was not to duplicate the published 
ion-chamber surveys. Instead, the objective was to take a few selected probe locations 
and obtain data for a wide range of propellant flow rates and ion-chamber configurations. 
Probe A (fig. 7) was selected as a representative sample of the primary-electron region. 
Probe B was located on the ’’critical field line” that separates the primary- electron 
region from the rest of the chamber. (The name "critical field line” was apparently 
first used by King, Poeschel, and Ward, ref. 119. ) Probe C was selected as a repre- 
sentative sample of the ion chamber outside of the primary- electron region. 

The ion-chamber probe data were obtained at several emission currents for each of 
several mass flows. Typical data from probe A are shown in figure 17 for one propel- 
lant mass flow. The density of Maxwellian electrons first increased with discharge, 
reached a maximum near the "knee” of the curve for ion-beam current, and then slowly 
decreased at higher discharge currents. The ratio of primary to Maxwellian electron 
densities increased continuously with discharge current throughout the range investigated. 
The electron temperature also increased with discharge current, but much less rapidly. 

The probe data were cross plotted to determine ion-chamber parameters for 1000- 
eV-per-ion condition defined as maximum propellant utilization. For the data in fig- 
ure 17 the value of 1000 eV per ion is reached at a discharge current of about 2. 1 am- 
peres. These cross-plotted values are presented in figure 18. The scatter in the data 
of figure 18 is to a large extent a measure of the uncertainties involved in obtaining and 
analyzing ion-chamber probe data. The overall uncertainty was estimated by Strick- 
faden and Geiler (ref. 21) as being about 30 percent, with most of the error associated 
with the straight-line assumption of primary- electron contribution. The frequent re- 
placement of probes was necessary because of cathode sputtering and no doubt added to 
the total scatter. This sputtering problem could have been avoided by using a hollow or 
oxide cathode. A hollow cathode would have introduced a perturbation in neutral density 
and thus analytical uncertainty. An oxide cathode would be likely to affect utilization 
measurements at very high discharge losses through evaporation from ion bombardment. 

The temperature of the Maxwellian electrons in the primary-electron region 
(probe A, fig. 18(a)) averaged about 5 eV, which agrees closely with values obtained 
previously at lower utilizations (ref. 123). This temperature drops slightly to about 
4 eV at the critical field line and is only about 2. 5 eV outside the primary-electron re- 
gion. 

The ratio of primary -to -Maxwellian electron densities averages about 0. 4 for the 
primary-electron region (probe A, fig. 18(b)). This ratio drops to about 0. 1 for the 
critical field line. Similar data are not shown for the ion chamber outside the primary- 
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Figure 17. - Typical variation of performance and probe data 
with discharge (anode) current. Probe A, 131 mA equiv. 
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(a) Maxwellian-electron temperature, probe A. 



(c) Maxwellian-electron temperature, probe B. 



(d) Primary-to-Maxwellian density ratio, probe B. 



(e) Maxwellian-electron temperature, probe C. 



(f) Probe B-to-A Maxwellian density ratio. 
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Propellant flow rate, mAequiv 


(g) Probe C-to-A Maxwellian density ratio. 

Figure 18. - Ion-chamber probe data for maximum 
utilization. Symbols same as figures 15 and 16. 
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electron region (probe C) because the ratio was too low to measure accurately (==0. 02). 

It should be noted that the primary -to-Maxwellian density ratio in the primary -electron 
region is well above the ^0. 1 value found previously at more normal operating condi- 
tions. This high value is no doubt due to the 1000-eV-per-ion discharge losses associ- 
ated with figure 18. One should also note that the electron parameters for this high dis- 
charge loss indicate that most of the ionization will take place in the primary-electron 
region, as assumed in the analysis. The mean energy of the primary electrons ranged 
from about 35 to 40 eV at probe A, ranged from 30 to 35 eV at probe B, and was less 
than 30 eV at probe C. 

The data shown in figure 18 show no significant trends with propellant mass flow or 

ion- chamber surface area. The use of the constant values of 5 eV for Maxwellian - 

electron temperature and 0. 4 for primary-to-Maxwellian density ratio appears justified 

for the primary-electron region values of equation (26). ^ These values, together with a 

value of 0. Oil meter for V /A (which can be obtained from a magnetic -field map such 

P P 18 

as shown in fig. 3), give a neutral density of 3. 11x10 per cubic meter. This neutral 
density, equation (4), an equivalent sharp-edged orifice area of 0. 00196 square meter, 
and an assumed wall temperature of 500 K result in a neutral loss rate equivalent to 
56 milliamperes. This value is quite close to the experimentally determined value of 
52 milliamperes (fig. 15). This close agreement is well within the possible error of the 
analytical model and constitutes an independent check on the analytical model. 


Comparison With SERT II Data 

It is of interest to compare the results of this investigation with published data. 
Probably the most documented thruster of divergent-field design is the one used on the 
SERT II flight (ref. 110). The SERT II throttling data (ref. 109) were reviewed and 
found to show the same nearly constant loss rate for un-ionized propellant that was found 
in this investigation. The comparison was not possible for the full throttling range above 
250 eV per ion. But this value was high enough to be in the "knee” of the loss curve and 
thus be substantially proportional to maximum utilization as defined in this report. The 
un-ionized propellant loss was equivalent to a current range from 45 to 50 milliamperes 
for 250 eV per ion and a total propellant flow range of 2. 5 to 1. (In a private communi- 
cation with Byers has stated that fig. 6(b) of ref. 109 is incorrect. The utilization at 
250 eV/ion and a beam current of 0. 20 was close to 0. 80. ) The level of un-ionized pro- 

^The temperature enters as a square root in eq. (26), so that the possible variation 
of temperature indicated by fig. 18(a) would enter as less than a linear effect. In a simi- 
lar manner, the possible variation of primary-to-Maxwellian density ratio would have 
less than a linear effect over the range of interest (see fig. 4). 
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pellant loss rate predicted from the theory of this investigation, however, was substan- 
tially above the observed equivalent range of 45 to 50 milliamperes. An unusual aspect 
of SERT II ion-chamber operation (with the accelerator system used for the throttling 
data, ref. 109) was for formation of several -volt ridges in plasma potential near the 
anode (refs. 85 and 113). Although the cause of these potential ridges is not clear, their 
presence probably reduced radial ion losses in their vicinity. If one neglects the bound- 
ary area of the primary-electron region near these ridges, the calculated neutral loss 
becomes much closer to the experimental values. Also, the introduction of propellant 
through the hollow cathode (about one-third of the total propellant flow at SERT II design 
conditions) was a departure from the uniform conditions assumed for the theory deriva- 
tion. Thus, although SERT II data are in good qualitative agreement with the theory of 
this investigation, it is apparently possible to obtain special plasma conditions that will 
affect the quantitative agreement of experimental data with the theory presented in this 
report. There is no evidence for reduced radial ion loss in the present investigation, 
inasmuch as probes A and B had an average potential difference of less than 1 volt. 


CONCLUDING REMARKS 

The most significant results of this investigation are the analyses of maximum pro- 
pellant utilization and the experimental agreement that was obtained with these analyses. 
The significance of these results is in the implied ability to predict maximum utilization 
of new designs from the field shape, and the inverse ability to select field shapes that 
should give high utilization. 

The results of this investigation can be expressed simply: the loss rate of un- 
ionized propellant at maximum utilization is nearly a constant over a wide range of total 
propellant flow rate. The simplicity of this result permits rapid prediction of total flow- 
rate effects from operation at one total flow rate. 

In addition to the general design application, the results of this investigation are 
particularly significant for throttled operation of a fixed -configuration thruster. Because 
of variations in available power, many missions require considerable operation at less 
than maximum thrust. Further, to avoid either large power losses or large propellant 
losses, such part-thrust operation should be near the "knee" of the ion-chamber loss 
curve. This knee will shift substantially in the same manner as maximum propellant 
utilization. Inasmuch as the un-ionized-propellant loss rate will be a larger fraction of 
total flow rate at part-throttle operation, one can expect the most difficult requirement 
for utilization to be at the lowest thrust level to be used. Without any configuration 
changes accompanying thrust reduction, then, a thruster with good part-thrust perform- 
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ance would be expected to have very high utilization at full thrust. Good qualitative 
agreement of these results was found with published SERT II throttling data. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 17, 1971, 

113-26. 
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